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Abstract. In this paper we compare data to theory. We 
use a compilation of the most recent cosmic microwave 
background (CMB) measurements to constrain Hubble's 
constant h, the baryon fraction fib, and the cosmologi- 
cal constant A Q . We fit h-, fib- and A D -dependent power 
spectra to the data. The models we consider are flat 
cold dark matter (CDM) dominated universes with flat 
(n s = 1) power spectra, thus the results obtained ap- 
ply only to these models. CMB observations can exclude 
more than half of the h — fib parameter space explored. 
The CMB data favor low values of Hubble's constant; 
h w 0.35. Low values of fib are preferred (fib ~ 0.03) 
but the x 2 minimum is shallow and we obtain fib < 0.28. 
A model with h w 0.40, fib ~ 0.15 and fl c dm ~ 0.85 is per- 
mitted by constraints from the CMB data, BBN, cluster 
baryon fractions and the shape parameter T derived from 
the mass density power spectra of galaxies and clusters. 

For flat-A models, the CMB data, combined with 
BBN constraints exclude most of the h — \ plane. Models 
with fl a as 0.3, A « 0.7 with h Rt 0.75 are fully consistent 
with the CMB data but are excluded by the strict new q Q 
limits from supernovae (Perlmutter et al. 1997). A com- 
bination of CMB data goodness-of-fit statistics, BBN and 
supernovae constraints in the h — X plane, limits Hubble's 
constant to the interval 0.23 < h < 0.72. 

Key words: cosmic microwave background — cosmology: 
observations 



1. Introduction 

A new technique is coming on-line producing a small rev- 
olution in our ability to evaluate cosmological models 
(Bond et al. 1994). Measurements of fluctuations in the 
cosmic microwave background (CMB) over a large range of 
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angular scales have become sensitive enough to distinguish 
one model from another. This technique is truly cosmolog- 
ical and independent of previous methods. It probes scales 
much larger and times much earlier (z » 1000) than more 
traditional techniques which rely on supernovae, galaxies, 
quasars and other low-redshift objects. 

Our ignorance of h, fib and A D is large and there are 
many inconsistent observational results. For example, af- 
ter 60 years of effort the favorite Hubble constants of re- 
spected cosmologists can differ by more than a factor of 
2 (0.40 Si h < 0.90) with error bars small compared to 
this interval. Thus, information on Hubble's constant and 
other cosmological parameters from the independent and 
very high redshift CMB data is important. Over the next 
decade observations at small angular scales have the po- 
tential to determine many important cosmological param- 
eters to the ~ 1% level (Jungman et al. 1996). We expect 
CMB measurements to tell us the ultimate fate of the Uni- 
verse (fl ), what the Universe is made of (fib, ^cdm) and 
the age and size of the Universe (h) with unprecedented 
precision. 

The COBE detection of temperature fluctuations in 
the CMB (Smoot et al. 1992, Bennett et al. 1996) has 
constrained the amplitude and slope of the power spec- 
trum at large angular scales. Scott et al. (1995), using a 
synthesis of the most current CMB measurements at that 
time, demonstated the existence of the predicted acoustic 
peak. Using compilations of CMB measurements, Ratra et 
al. (1995), Ganga et al. (1996), Gorski et al. (1996), White 
et al. (1996) have compared groups of favored models to 
the data and obtained interesting constraints on fl , n and 
the normalizations of various flat and open CDM models. 

The most recent ground-based and balloon-borne ex- 
periments (Netterfield et al. 1997, Scott et al. 1996, Piatt 
et al. 1996, Tanaka et al. 1996) are providing increasingly 
accurate CMB fluctuation measurements on small angular 
scales. In this work we take advantage of these new mea- 
surements and of a fast Boltzmann code (Seljak & Zal- 
darriaga 1996) to make a detailed exploration of two large 
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regions of parameter space: the h — Qb and h — \ planes. 
Such an approach allows us to place independent CMB- 
derived constraints on h, Qb and A based on goodness-of- 
fit statistics. 

Our results are valid only in the context of the models 
we have considered; we assume COBE normalized, Gaus- 
sian, adiabatic initial conditions in spatially flat universes 
(k = 0, Q + A = 1) with Harrison-Zel'dovich (n s = 1) 
temperature fluctuations. For Gaussian fluctuations the 
power spectra, Ce, uniquely specify the models. Our lim- 
its include estimates of the uncertainty due to the COBE 
normalization uncertainty as well as the Saskatoon abso- 
lute calibration uncertainty. We use h = H o /10Q km s _1 
Mpc -1 . 

Our CMB-derived constraints on h, Qb and A Q are in- 
dependent of BBN and other cosmological tests. Where 
the CMB constraints are not as tight as other methods, 
the existence of overlapping regions of allowed parameter 
space is an important consistency check for both; unknown 
systematic errors can be uncovered by such a comparison. 

In Section 2 we describe the x 2 calculation. In Section 
3 we provide an overview of the physics of fluctuations and 
the power spectrum models used in the fit. In Sections 4 
and 5 we discuss our results for the h — fib and h — X Q tma* 
planes respectively and combine them with a variety of I(Wt) — 2tt 
other cosmological constraints. In Section 6 we summarize ^=2 
and discuss our results. 



a a (N). The models are indexed in the 2-D parameter 



space by i and j. The band power estimates of the data 
and models are defined respectively as 
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where ST°^ S (N) is the rms temperature fluctuation ob- 
served by the Nth experiment, We(N) are the experiment- 
specific window functions (White & Srednicki 1995) and 
the deconvolving factors I(We(N)) are the logarithmic in- 
tegrals of the We (Bond 1995) defined as 
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2. Data Analysis 

2.1. The Recipe 

• assemble all available CMB experimental results with 
the corresponding window functions 

• determine the region of parameter space one would like 
to explore (we choose the h — Oj, and h - \ a planes for 
spatially flat models) 

• use a fast Boltzmann code (Seljak & Zaldarriaga 1996) to 
obtain the power spectra for a matrix of models covering 
the desired region of parameter space 

• convolve each model with the experimental window 
functions and fit the result to the data by producing a 
X 2 surface over the chosen parameter space 

• compare the results with other cosmological constraints 

2.2. The Equations 

We calculate the x 2 surface over a matrix of models 



where £ max = 1200 since for all N, W e (N) w for t > 
1200. The angular scale probed by the Nth experiment is 
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and use it to indicate the regions of parameter space pre- 
ferred by the data. The sum is over the CMB detections 
plotted in Figure |l| and listed in Table 1 in the form of 
flat band power: STf ata (N) with the corresponding errors 



In equation (2), the convolution of the unknown real 
power spectrum of the CMB sky, Cereal), with Wi is 
the observed temperature rms, ST°^ S . One cannot com- 
pare the rms results from different experiments with each 
other unless the influence of the window function has been 
removed, i.e., deconvolved. The division by I(We) is this 
deconvolution and works optimally when £(£ + l)Ce(real) 
is a constant across the range of £ sampled by the experi- 
ment. Equation (4) has been written suggestively to clar- 
ify this deconvolution. Notice that the models are treated 
like the real sky: first they are convolved with the window 
function and then the division by I (We) deconvolves the 
window function. Thus our model points are "convolved- 
deconvolved" power spectra. Setting Ce = a/£(£ + 1) in 
equation (2) or (3) yields ST lcff = (a/27:) 1 / 2 = [£(£ + 
l)Ce/2n] 1 ^ 2 which is the origin of the units of the y-axis of 
Figure |l| and why it is reasonable to plot the input model 
and the flat band power estimates on the same plot. 

Using the data and equation (3) in equation (1) yields 
a x 2 value for every point in the matrix of models. Figure 
is a picture of the ingredients of our x 2 calculation. The 
CMB measurements from Table 1 are plotted. The thick 
solid line is an h — 0.30 model with fib = 0.05. The large 
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Fig. 1. The power spectrum of observed CMB temperature 
fluctuations (from Table 1) as a function of angular scale. The 
FWHM of the window functions are indicated by the horizon- 
tal lines at each data point. Two representative models are also 
plotted; both have fit = 0.05. The differences between the con- 
volved h — 0.3 model (large circles) and the measurements is 
used to calculate a single point of the \ 2 surface shown in Fig- 
ure [s| The small boxes above and below the Saskatoon points 
indicate the ±14% absolute calibration uncertainty shared by 
these 5 points. The dotted line is a 5th order polynomial fit to 
the data. On the upper axis, Ofwhm = (32ln2) 1 ^ 2 sin~ 1 ( 2 ^ +1 ) 
(Bond 1995). 



open circles are the convolution of this h = 0.30 model 
with the experimental window functions. The difference 
between the resulting "convolved-deconvolved" points and 
the measurements is used to calculate the x 2 values. By 
comparing the data with these "convolved-deconvolved" 
points rather than with the direct [£(£+l)Ce/2n] 1 / 2 values 
of the original h = 0.30 model, we are accounting for the 
experimental window function even in the regions of sharp 
peaks and valleys. The discrepancy between the open cir- 
cles and the original h = 0.30 model (solid line) is a mea- 
sure of the necessity of this convolution-deconvolution pro- 
cedure; it is unnecessary except in the sharp peaks and 
valleys. The dotted line is a fifth order polynomial fit to 
the data points and is used in Figures ^| and to repre- 
sent the data. Notice that the amplitude of the primary 
acoustic peak in the data is A pea k ~ 80 /iK. 

Most experimental results contain slightly asymmetric 
error bars and, most commonly, the upper error bar is 



Fig. 2. Representative power spectra showing the h- and 
fib-dependence. In each panel h is fixed while ilb takes on 
the values indicated. The largest values of Qj, have the largest 
Doppler peaks. Notice that as h increases, A pea k increases for 
large Qb but decreases for small Qb', thus at high h the peak 
height is an excellent baryometer. The dotted line is the poly- 
nomial fit to the data in Figure |l| and is the same in all pan- 
els. All models are spatially flat (k = 0, A D = 0), Harrison 
Zel'dovich (n 3 = 1) normalized to the COBE 4-year results at 
1 = 6. 



larger. For a given data point, this asymmetry constrains 
models lower than the data point a bit more rigidly than 
it does models which are higher than the data point. To 
approximate this asymmetry for the \ 2 calculation, we 
toggle the error bar depending on whether the model is 
above or below the data point. 



3. CMB Power Spectra 

3.1. General Features 

Figures ^| and [5| are samples from the matrices of power 
spectra covering the parameter space explored. They show 
the influence of h, flh and A Q on the angular power spec- 
trum of CMB fluctuations. Generic features are: 

• a flat Sachs- Wolfe plateau for £ <> 30 due to the grav- 
itational potentials of superhorizon-sized density fluctua- 
tions. 

• a primary acoustic peak at £ pea k ~ 250 accompanied by 
secondary peaks and valleys for £ J; 2 £ pea k 



4 



Charles H. Lineweaver et al.: Constraints on h, fib and A from Cosmic Microwave Background Observations 



• a cut-off at I <; 1000 due to averaging along the line 
of sight as the photons traverse the finite thickness of the 
surface of last scattering 

The physics of the primary acoustic peak is the most 
relevant for our purposes since this peak dominates the 
fits. At sub-degree angular scales acoustic oscillations of 
the baryon-photon fluid at recombination produce peaks 
and valleys in the CMB power spectrum. We may un- 
derstand their general features by employing the driven 
oscillator interpretation of Hu (1995) and Hu & Sugiyama 
(1995a). The waves are set up by the pressure of the pho- 
tons, which are tightly coupled to the baryons and subject 
to the driving force of the gravitational potential. Ignoring 
the neutrinos, there are three particle species present at 
recombination (baryons, photons and dark matter) which 
determine two fundamental ratios: the baryon-to-photon 
ratio, varying as flbti 2 , and the dark matter-to-photon 
ratio, varying as fl c d m h 2 . We expect the physics at this 
epoch to be largely invariant to changes which leave these 
quantities unaltered, i.e., we expect flbh 2 and fl c d m h 2 to 
control the shape of these oscillations and in particular 
the height and location of the primary peak A pea k , l pea k ■ 
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Fig. 3. Contours of A pea k \jj.K]. Although no data were used 
to make this plot, the 80 pK contour is thick to indicate the 
approximate peak height of the data. For large h, the strong 
^-dependence of the peak heights is responsible for the nar- 
row spacing between contours and for the narrow range of fit 
permitted by the data (see Figure []). As h decreases, A pea k 
loses its S7b-dependence thus permitting a wider range of fit,. 
The grey band is the region favored by big bang nucleosyn- 
thesis: 0.010 < fl b h 2 < 0.026 (e.g. Copi et al. 1995, Tytler & 
Buries 1997) . All models in this plot have fi D = fib + flcdm = 1 
and A = 0. 
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Fig. 4. Contours of £ pea k for the same models as in Figure |3[ 
The thick 280 contour is indicative of the £ pea k favored by the 
data. In the lower right, £ pea k is too small to suit the data. 

3.2. Effects ofil b h 2 and fl c d m h 2 
3.2.1. Peak Height 

Figures 0, ||, || are a triptych; they show three different 
ways of looking at the h — fib models we have explored. 
Figure shows the entire power spectra but only for a 
sparse sample of parameter space. Each panel shows 5 
samples from a vertical strip of Figures | and |. Figures 
U and ^ fully sample the matrix of models but with only 
one-parameter characterizations of the spectra. 

It is important to distinguish aspects of the x 2 surface 
produced by the underlying models from the aspects pro- 
duced by the data. The data currently available are in the 
range 2 t e ft 600 and are insensitive to most of the 
fine details in the models seen at I <; 2 £ pea k- Since the 
Doppler peak of the spectra is the most prominent feature, 
its amplitude A pea ^ is the simplest 1-parameter character- 
ization of the the relative topology of the parameter space 
and is an excellent tracer of what the yf 2 contours will look 
like. 

Figure || is a contour map of A pea k created by plotting 
the maximum amplitude of each model power spectrum. 
Although no data are involved, the A pea k — 80 /iK contour 
is thick to indicate the approximate peak height of the 
data in Figure [l]. If A pea k were the only important feature 
relevant for the fit, the \ 2 minimum would straddle the 
A P eak ~ 80 fjK contour and be contained by boundaries 
parallel to these peak contours. This is almost the case 
(see Figure ||) . 

In Figure || one can see that for a given h, A pea k in- 
creases monotonically with f2{,: A pea k is low at the bottom 
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of the plot and high at the top. The grey band marks the 
BBN region bounded by contours of £l&/i 2 = constant. 
In the upper right, higher values of fibft 2 lead to larger 
Doppler peaks due to the enhanced compression caused 
by a larger effective mass (more baryons per photon) of 
the oscillating fluid. 

If we move along the BBN region from right to left, 
Apeak increases. Since Q{,h 2 — constant in this region, 
we are seeing the influence of n, c d m h 2 . As h goes down, 
Qcdmh 2 goes down, and therefore A pea k goes up. This 
can be understood more physically. Lowering h increases 
the fraction of the total energy density locked up in the 
photons. As this component cannot grow in amplitude (it 
is oscillating), it retards the potential evolution with re- 
spect to purely matter-dominated fluid - in other words, 
it causes a decay in the gravitational potential. An oscilla- 
tor driven by a decaying force will actually obtain larger 
amplitudes around its zero point than one subject to a 
constant force, assuming the same initial conditions. At 
fixed flbh 2 , this leads to a larger primary acoustic peak. 

Contours parallel to the BBN region with A pea k in- 
creasing to the upper right, are imposed by f^/i 2 . Nearly 
vertical contours are imposed by Q c dmh 2 {£l c dm ~ f2 D = 1) 
with Ap ea k increasing to the left. In summary, Q, h 2 and 
Qcdmh 2 produce competing effects; when Q h 2 ft A pea k ft 
and when £l c d m h 2 ft A pea k -IJ-. The two effects combined 
explain the shape of the contours in Figure |[ 

In Figures [| and || we see the effect of varying A and 
h with fl h 2 held fixed. Since fl h 2 is fixed, for any given 
A Q = constant, when h ft Q c d m h 2 ft and thus A pea k ft. 
Also since fl + Sl c d m + A„ = 1, for any given h when A ft 
we have £l cd mh 2 ft thus A pea k ft. We can see both of these 
effects in the A pea k contours of Figure ||. As in Figure || no 
data were used to make this plot and the 80 fiK contour is 
thick to indicate the approximate peak height of the data. 
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Fig. 5. Representative power spectra showing the h- and 
Ao-dependence. In each panel h is fixed while A takes on 
the values indicated. The largest values of A have the largest 
Doppler peaks. This is predominantly the effect of fi c d m /i 2 
decreasing since the baryon-to-photon ratio has been fixed: 
fit/i 2 — 0.015 and (f2 c dm + fit + A = 1). Notice that as h 
increases, A pea k decreases and £ pea k shifts to larger scales. The 
dotted line is the polynomial fit to the data in Figure |l| and 
is the same in all panels. All models are spatially flat (k = 0), 
Harrison Zel'dovich (n s — 1) and normalized to the COBE 
4- year results at £ = 6. 



3.2.2. Peak Position 

Hu & Sugiyama (1995b) identify three factors which, when 
increased, increase £ pea k- Qbh 2 , A and (1 — fl — A D ). And 
one factor which, when increased, decreases £ pea k, Q h 2 . 
For all our models 1 — £1 — A Q = 0, so this focusing effect 
plays no role here. The effect of Vlbh 2 can be understood 
using its relationship to the sound speed at recombination: 
c s = (1 + af^ft 2 ) -1 / 2 where a is a constant, thus if fibft 2 ft 
c s ft. The size of the fluctuations is controlled by the sound 
horizon at decoupling, thus 9 pea k oa c s . If c s ft then pea k ft 
and thus tpeak ft- We can see this clearly in the upper right 
of Figure [|. Qbh 2 and fl c dmh 2 have opposite effects on the 
peak amplitude and location. When Vtbh 2 ft l pea k ft and 
A P eak ft- When VL cdm h 2 ft f peak ft and A peak ft. 

The location of the peak plays an important role in pre- 
ferring the low h region of iso-A pea k contours. The CMB 
data have £ pea k ~ 280 which corresponds to an angular 
scale of ~ 0°.5. In the lower right of Figures [| and (?| 



ipeak ~ 210 (substantially lower than 280). Thus the data 
disfavor these high h models. 

There is an interesting apparent inconsistency in Fig- 
ure [?]. In the entire figure f^ft 2 = 0.015 = constant. Along 
the dashed line fl h 2 = n b h 2 +fl cdm h 2 = 0.015+(0.41) 2 = 
constant. Thus as we follow the dashed line up, A Q in- 
creases and we should see £ pea k increasing since, as we 
mentioned earlier if A ft then £ pea k ft- This is not the case 
however; as we follow the dashed line upwards, £ pea k de- 
creases. A scaling argument can clarify this. The angle 9 
subtended by an object of size ct at an angular distance 
dang is 6 ~ ct/d ang . Thus, 

t k 1 d ang (h,n ,\ ) h- l f{n ,\ ) /(A n ) | 

9peak C tree (^o^- 2 ) ^ 2 ^ 

where 9 peak is the angular scale of the Doppler peak. The 
physical scale of the peak oscillations is some fixed fraction 
of the horizon oc c t rec . The time of recombination scales 
as (floh 2 )^ 1 / 2 (which is constant along the dashed line) 
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and the angular distance d ang = h 1 f(fl , X ). Since in 
our flat models tt = 1 — A , f(Sl , A D ) — > /(A c ) where 



/(Ac) 



dz 



\(l + zf-\ {l + z)- 



Hi/2" 



(7) 



which, when inserted into equation (6) is an expression 
of the monotonic relation, when A Q f|~ then £ pea k fl- How- 
ever, for flat A models with fl^h 2 and £l c dmh 2 fixed, one 
cannot change A without changing h. So, in the case 
we are considering equation (6) is telling us that when 

f ( A ) 

1 y h °' -]\ £p e ak 1\- The h scaling can be understood as the 
effect of larger universes: a given physical size at a larger 
distance subtends a smaller angle. 

Following the pioneering work of Hu & Sugiyama 
(1995a, 1995b), in this section we have presented contours 
of the solutions of the Boltzmann equation and we have 
discussed how the behavior of A pea k and £ pea k can be ex- 
plained in terms of Q,bh 2 and fl c d m h 2 . These contour plots 
are particularly relevant for the next section where we de- 
scribe the regions of solution space preferred by the CMB 
data. 



Fig. 6. Contours of A peak [fiK] with Q b h 2 = 0.015. As in Fig- 
ure^ no data were used to make this plot. The 80/iK contour is 
thick to indicate the approximate peak height of the data. For 
a given h, as Ao increases, A pea k increases. The grey area marks 
the permitted area (la) from new supernovae results (Perlmut- 
ter et ai. 1997). All models in this plot have Q + A = 1. 
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Fig. 7. Contours of £ pea k for the same models as in Figure ^. 
Although no data were used to make this plot, the thick 280 
contour is in the center of the region favored by the data. The 
dashed line is discussed in Section 3.2.2. 




4. h — fit, Results and Discussion 

The CMB data can constrain any parameter that changes 
the power spectrum at a level comparable with the error 
bars on the data. Figure || displays the \ 2 contours for 
models in the h — Clb plane and contains one of the main 
results of this paper. CMB observations alone exclude at 
> 95% CL, more than half of this parameter space. The 
solid and dotted contours in Figure || are goodness-of-fit 
contours labeled with the probabilities of finding a \ 2 less 
than the calculated value at that point. For example, in 
our case there are 21 degrees of freedom (24 data points - 2 
fitted parameters - 1 COBE normalization). The probabil- 
ity of obtaining a x 2 value less than 23.8 is 68.3 % assuming 
uncorrelated measurements and Gaussian errors. Thus we 
have labeled the \ 2 = 23.8 contour '68'. The levels plotted 
are x 2 = [17.2,23.8,32.6] corresponding to 30%, 68% and 
95% respectively. 

The solid-line contours include our estimate of the 14% 
absolute calibration uncertainty shared by the 5 Saskatoon 
points. The thick dotted (68%) and thin dotted (95%) 
contours do not include the Saskatoon calibration uncer- 
tainty. Thus the Saskatoon errors do not affect the 95% 
contour in the left and upper right, nor the left side of 
the 68% contour. The effect of the Saskatoon calibration 



uncertainties is discussed further in Section 4.1 



Any data sensitive only to the peak height will yield 
contours as seen in Figure pi Figures | and | are similar 
so peak height is dominating the fit. If A pea k were the 
only factor, then a region straddling the A pea k ~ 80 
contour would be equally preferred by the data. However 
the goodness-of-fit statistic prefers a region that does not 
follow exactly any iso- A pea k contour. There is a preference 
for the low h part of an iso-A pea k region. 
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Lower h is preferred because of the location of i pea k 
in the data. One can see in Figure that at smaller h 
the position of the peak is more aligned with the dotted 
(data) line. This can also be seen in Figure ^ where the 
(■peak values in the lower right are substantially less than 
the ~ 280 of the data. This £ pea k mismatch disfavors the 
high h part of the preferred iso-^4 pea fe region. 




.20 .30 .40 .50 .60 .70 .80 .90 1.0 
h 



Fig. 8. Likelihood region and goodness-of-fit contours of y 2 . 
The similarity of this figure and Figure ^ is a measure of the 
domination of the fits by A pea k- The solid and dotted lines 
are goodness-of-fit contours. The solid contours (30%, 68%, 
95%) include our estimate of the Saskatoon calibration errors; 
the dotted contours (68%(thick), 95%(thin)) do not. The lower 
thin 95% dotted contour is obvious but the upper 95% dotted 
contour is identical to, and therefore hidden by, the solid 95% 
contour. The two dark grey areas are the 68% likelihood regions 
for SkO (left) and Sk— 14 (right) (see text). Small values of h are 
preferred. The region preferred by big bang nucleosynthesis 
(0.010 < tt b h 2 < 0.026) is the light grey band. The small box 
and triangle at Qj, w 0.05 are the two model spectra plotted in 
Figure |l[ 



4-1- Saskatoon Calibration Error 

The 5 Saskatoon (Sk) measurements (Netterfield et 
al 1997) apparently span the Doppler peak and play an 
important role in our fitting procedure. All 5 points share 
a 14% absolute calibration error as indicated by the small 
squares above and below their central values in Figure [l| 
There is no systematic way to handle systematic errors. 
Adding 14% errors in quadrature to the statistical errors is 



not appropriate because these errors are 100% correlated 
in the 5 Sk points. We have examined the effect of these 
correlated errors is several ways. We have made x 2 surfaces 
from the CMB data: 

• using the 5 Sk points as listed in Table 1 ( "SkO" ) 

• adding 14% to the 5 Sk points, i.e., using the high small 
squares in Figure |l| as the central values ("Sk+14") 

• subtracting 14% from the 5 Sk points, i.e., using the 
lower small squares in Figure [l] as the central values 
("Sk-14"). 

The minimum \ 2 f° r these three cases are respectively, 
20, 33 and 13 corresponding to goodness-of-fit contours of 
48%, 96% and 9%. Thus SkO gives a reasonable fit, Sk+14 
gives a bad fit and Sk— 14 gives a fit that is a bit too good. 

We take a conservative approach in our treatment of 
these calibration errors. We create a new surface by adopt- 
ing at each point in the plane the minimum of these three 
X 2 surfaces (Sk+14, SkO and Sk-14). In Figure | we plot 
the iso-% 2 contours from this new surface. We will call 
this the "Sk+14" case. For this new surface the Sk— 14 
X 2 surface determines the contours in the lower right while 
the SkO x 2 surface determines the contours to the left 
and upper right. This is easy to understand with Figure 
||; SkO has a higher peak than Sk— 14. Thus SkO prefers 
Apeak ~ 80/LtK and Sk-14 prefers A peak ~ 70 fjK. The 
Sk+14 surface is effectively eliminated because of its poor 
goodness-of-fit. 

The contours in Figure ^ are conservative in the sense 
that they encompass a larger region than the correspond- 
ing contours without including the Saskatoon calibration 
uncertainty. The dotted contours are the 68% and 95% 
goodness-of-fit levels from SkO. The solid lines are the 
goodness-of-fit levels from our Sk±14% method. Thus Fig- 
ure |^ shows where and by how much the Saskatoon cali- 
bration uncertainty reduces the constraining ability of the 
CMB data. 

It is interesting to note that the rest of the data thinks 
that the Saskatoon values are too high in the sense that the 
X 2 surface from Sk— 14 has a lower minimum than the SkO 
surface. This can be quantified; we can guess at the correct 
Sk calibration by using the non-Sk data and the models 
while treating the common calibration uncertainty of the 
Sk points as a free parameter, i.e., we let the calibration 
vary and assume that one of the models is correct. The 
lowest x 2 is obtained if the Saskatoon points are reduced 
by 24%; the non-Sk data would like to see the Sk points 
all lowered by 24%. This preference is however shallow: 
X 2 min (Sk-2A%) = 11.4 while X 2 min (Sk-U%) = 12.7 and 
XL„(^-0%) = 19.9. 

The fact that the fit to Sk— 14 is a bit too good can be 
interpreted as some combination of the following: 

• just chance 

• non-Sk data prefer a peak height lower than the Sk mea- 
surements 

• the error bars of the non-Sk points are over-estimated 

• detections of small significance do not get reported or 
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Fig. 9. Same as Figure g except instead of the Saskatoon un- 
certainty, here we show the effect of the COBE normalization 
uncertainty. The dotted contours are the same as in Figure 
|§| The solid contours include our estimate of the effect of the 
COBE normalization uncertainty. They are the contours of a 
X 2 surface constructed from the minima of three x 2 surfaces, 
one each for Q = 16.4, 18.0 and 19.6. The higher Q value lets 
Qb be as high as 0.28. The limits on h do not change substan- 
tially and thus small values of h are still preferred independent 
of the COBE normalization uncertainty. 

are reported as upper limits which have not been included 
in the data set 

• observers are finding what they are supposed to find. 

Reasonable goodness-of-fit is a prerequisite for 
X 2 minimum parameters that are meaningful. In addi- 
tion to the goodness-of-fit contours presented thus far 
we have used the minima of the SkO, Sk— 14 and Sk+14 
X 2 surfaces to define the best fit parameters and to define 
confidence levels around these minima. This procedure is 
correct when the goodness-of-fit is within some plausible 
range, commonly [5%, 95%], thus the Sk+14 minimum is 
of questionable value. To obtain likelihood intervals on pa- 
rameter values at the minimum of a 2-D \ 2 surface one 
takes Xmin + 2-3 as the 68% contour and Xmin + 6-2 as 
the 95% contour (e.g. Press et al. 1992). In Figure |, the 
two dark grey areas are the 68% likelihood regions for SkO 
(left) and Sk— 14 (right). The preferred h value is low and 
the confidence levels are smaller than the corresponding 
goodness-of-fit contours. 

4-2. Normalization Uncertainty 

So far we have assumed the COBE normalization 
Q\n s =i — 18 fiK at £ = 6. However there is a 1.6 /zK 



overall uncertainty on this value. To include this uncer- 
tainty in our results we make x 2 contours for Q = 16.4, 
18.0 and 19.6 fiK. The minimum of the surfaces (produced 
by a procedure analogous to that used to make Figure ||) 
is shown in Figure ^ where the dotted contours do not, 
and the solid contours do, include the COBE normaliza- 
tion uncertainty. The normalization uncertainty increases 
the upper limit on Q, b from 0.17 to 0.28 (68% CL). The 
limits on h are robust to the normalization uncertainty. 
The lower right is still excluded. 

4-3. Sensitivity to "Outliers" 

One possible danger in using x 2 goodness-of-fit contours 
is the influence of outliers. The overall level (goodness- 
of-fit) and possibly the shape of the contours around the 
minimum can be controlled by outliers. This depends on 
whether the models differ much at the £ e ff of the sus- 
pected outlier. For example in Figure |, the COBE point 
at I — 25 is an "outlier" which does nothing more than 
raise the entire x 2 surface, since all models are the same at 
I ~ 25. This is not necessarily the case for the low MAX 
point (MAX HR) and certainly not for the 5 Saskatoon 
points. When the MAX HR point is excluded from the 
SkO surface calculation, the goodness-of-fit of the mini- 
mum improves from Xmin = 20 to Xmin = 15. The size 
and shape of the resulting contours do not change signif- 
icantly, probably because at £ w 140 the models are not 
substantially different. 




Fig. 10. Same as Figure g except we have added two more 
constraints: Clusters: 0.04 < Q b h 3/2 < 0.10 (White et al. 1993) 
and the galaxy and cluster power spectrum shape parameter: 
0.20 < T < 0.30 (Peacock & Dodds 1994). 
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4-4- Result Summary 

Our main results for the h — Of, plane are summarized 
in Figure ||. In the context of flat CDM universes, CMB 
observations can exclude more than half of the h — flb pa- 
rameter space explored. The CMB data favor low values 
of Hubble's constant; h ~ 0.35. A higher value such as 
h = 0.75 is permitted by the goodness-of-fit contours if 
0.05 < ilb < 0.09 but these values are not within the 1-er 
likelihood region and are excluded by big bang nucleosyn- 
thesis (BBN). Low values of fib are preferred (fib ~ 0.03) 
but the x 2 minimum is shallow and we obtain fib < 0.28 
(68% CL). 

These limits include our estimates for the Sk abso- 
lute calibration errors and the COBE normalization un- 
certainty. They are conservative in the sense that they are 
larger than the corresponding error bars around the pa- 
rameters at the minimum of the % 2 surfaces. Outliers do 
not seem to be biasing these results. 

The CMB and BBN are independent observations. 
Thus, the overlapping of regions preferred by BBN and 
the CMB did not necessarily have to exist. The fact that 
they do is a consistency test for these two pillars of the 
big bang scenario. 

The uncertainties in the BBN limits are dominated by 
systematic errors. Therefore we have used 'weak' BBN lim- 
its to try to avoid over-constraining the parameters. We 
have increased the limits of Copi et al. (1995), 0.010 < 
flbh 2 < 0.020, to include the higher values obtained by 
Tytler & Buries (1997), 0.022 < fl b h 2 < 0.026. The cred- 
ibility of the BBN constraints used here is important be- 
cause their incompatibility with the CMB constraints at 
h <; 70 is what excludes these models. If one takes the 
maximum upper limit on flbh 2 as 0.045 (Reeves 1994) then 
at the 68% CL all values of h between 0.30 and 1.0 are 
permitted by even the combination of the CMB and BBN 
limits. 

Figure [l0]is the same as Figure || except we have added 
two more constraints: the White et al. (1993) limits from 
cluster baryon fractions: 0.04 < flbh 3 / 2 < 0.10, and the 
CDM shape parameter from a synthesis of the power spec- 
tra of galaxies and clusters: 0.20 < T < 0.30 (Peacock & 
Dodds 1994). The similarity of the T limits and the SkO 
contours is interesting. The inconsistency of BBN limits 
with the White et al. (1993) limits (assuming fl = 1) is 
sometimes invoked as an argument for fl < 1. However, 
as Figure 10 shows, this inconsistency disappears for the 
low h values favored by the CMB data. An fl Q = 1 model 
with h w 0.40, fib ~ 0.15 and fl c dm ~ 0.85 is permitted 
by constraints from the CMB data, BBN, cluster baryon 
fractions and T. 



5. h — A Q Results and Discussion 

In the previous section we described results in the h — fib 
plane (k = 0, fl = 1, A G = 0). In this section we present 
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Fig. 11. Goodness-of-fit contours of \ 2 with flbh 2 — 0.015. As 
in Figure || the solid contours include the Saskatoon calibra- 
tion uncertainty, the dotted do not. The grey area marks the 
permitted area (1 o) from new supernovae results (Perlmutter 
et al. 1997). For all models in this plot, fl„ + A D = 1. 



results for the h — X D plane, (k = 0, fl D + A = 1). Figure 
pT| presents the x 2 contours from fitting the data to these 
models. We are exploring a plane orthogonal to Figure |[ 
The same procedure used to make the contours in Figure 
H was used here. The notation is also the same; the solid 
contours include the Saskatoon calibration uncertainties, 
the dotted do not. 

Figures |] and [ll] are similar for the same reason that 
Figures || and || are similar: A pea k is the dominant feature 
of the fit. There is a slight preference for the low h side 
of the preferred iso-A pea k contour. Figure |?] helps explain 
this preference. For a given A OJ i pea k ~ 280 is at a lower 
h than the Ap^ ~ 80 contour. 

In this plane we obtain 0.23 < h < 0.72 where 
these limits include our estimates for the uncertainties 
from the Saskatoon calibration, the COBE normalization 
and the BBN interval 0.010 < fl b h 2 < 0.026. Assum- 
ing fl b h 2 = 0.015, the CMB data yield A D <, 0.9 (68% 
CL) with lower values preferred. This is weaker than the 
traditional dynamical limit A c ^ 0.8. A standard flat- 
A model is fl ~ 0.3 and A ~ 0.7 with h w .70 This 
model is fully consistent with the CMB data. However 
the new SNIa results (Perlmutter et al. 1997) rule it out; 
-0.28 <, A <, 0.34 (la). 
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Fig. 12. Same as Figure [Ll| except we have added two more 
constraints and left out the dotted contours for legibility. Con- 
straints on the age of the universe (13.7 < t < 17.9) from the 
oldest stars (Bolte & Hogan 1995) are indicated by the thin 
solid lines. The dark grey area marks limits on the shape pa- 
rameter: 0.20 < T < 0.30 (Peacock & Dodds 1994). The thick 
dashed line through the minimum is fi cdm h 2 = (o.4i) 2 . 

6. Summary and Discussion 

We have explored solutions to the Boltzmann equation in 
two 2-D parameter spaces within the context of COBE 
normalized, Gaussian, adiabatic initial conditions in spa- 
tially flat universes with Harrison-ZePdovich temperature 
fluctuations. We have presented the topology which con- 
trols the shapes of the acceptable regions for the cosmolog- 
ical parameters h, fit,, and A Q . We have fit a compilation 
of the most recent CMB measurements to the models and 
have identified regions favored and excluded by the data, 
i.e., we have fit the Boltzmann solutions to the CMB data 
and obtained constraints on the cosmological parameters 
h, fib, and A Q . 

Figures ^| and |ll| contain the main results of this paper 
and are plots of the goodness-of-fit contours of the \ 2 from 
equation (1) using the data in Table 1. The solid contours 
mark the regions preferred by the data and include the 
Saskatoon calibration uncertainty. We obtain the follow- 
ing results: 

• Recent CMB data are precise enough to prefer distinct 
regions of parameter space and rule out most of the h — Qj, 
and h — A D planes at the 95% CL. 

• The CMB data favor low values of Hubble's constant; 
h ps 0.35. A higher value such as h — 0.75 is permitted 



by the goodness-of-fit contours if 0.05 < Q& < 0.09 but 
these values are not within the 1-a likelihood region and 
are excluded by big bang nucleosynthesis (BBN). 

• Low values of fib are preferred (fib ~ 0.03) but the 
X 2 minimum is shallow and we obtain fib < 0.28. 

• The CMB regions overlap with BBN. The fact that they 
do is a consistency test for these two pillars of the big bang 
scenario. 

• A fl a = 1 model with h « 0.40, fl b « 0.15 and 
flcdm ~ 0.85 is permitted by constraints from the CMB 
data, BBN, cluster baryon fractions and the shape pa- 
rameter r derived from the mass density power spectra of 
galaxies and clusters. 

• For flat-A models, the CMB data, combined with BBN 
constraints exclude most of the h — A Q plane. 

• Models with fl a w 0.3, A Q w 0.7 with h « 0.75 are fully 
consistent with the CMB data but are excluded by the 
strict new q Q limits from supernovae. 

• A combination of CMB data goodness-of-fit statistics, 
BBN and supernovae constraints in the h — \ Q plane, limits 
Hubble's constant to the interval 0.23 < h < 0.72. 

In this analysis we have made estimates of the errors 
associated with the Saskatoon absolute calibration uncer- 
tainty, the COBE normalization uncertainty, the uncer- 
tainty in BBN and we have considered the influence of 
possible outliers. We have assumed Gaussian adiabatic ini- 
tial conditions. We have conditioned on the values n s — 1 
and Yffe = 0.24. We have ignored the possiblity of open 
universes, tilted spectra, early reionization and any grav- 
ity wave contribution to the spectra. 

Our h result is 'low' and inconsistent with several other 
recent h measurements (0.65 ^ h 0.80) in the sense that 
our 68% goodness-of-fit contours do not, but our 95% con- 
tours do include these higher h values. The theoretical ad- 
vantages of a low Hubble constant have been presented in 
Bartlett et al. (1995). For example, if fl = 1 then a Hub- 
ble constant of 0.70 implies an age of 9.3 Gyr, younger 
than the estimated ages of many globular clusters. How- 
ever h ~ 0.40 yields an age of the universe of 16 Gyr, 
much more in accord with globular cluster ages. 

If new data show that A pea k ~ 70 rather than 
80 /iK and £ pea k ~ 230 rather than 280, then the range 
0.65 < h < 0.80 may be acceptable to both the CMB and 
BBN constraints. 
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Table 1. Data Used in the \ 2 Fit" and Plotted in Figure [j] 



Experiment 


reference 


4// 


° T *eff ±<J 


DMR1 


Hinshaw et al. (1996) 


3 


27.9^ 


DMR2 


Hinshaw et al. (1996) 


7 


24.61^ 


DMR3 


Hinshaw et al. (1996) 


14 




DMR4 


Hinshaw et al. (1996) 


25 


i r,+ 14.'4 

L -^-1.2 


FIRS 


Ganga et al. (1994) 


10 


29.4+™ 


Tenerife 


Hancock et al. (1996) 


19 


34.i±i2:H 


SP91 


Gunderson et al. (1995) 


60 


30.2t 8 5 j 


SP94 


Gunderson et al. (1995) 


60 


36.3+ji 3 ! 6 


Pythl 


Piatt et al. (1996) 


91 


54.0t&° 


Pyth2 


Piatt et al. (1996) 


176 


58.o±i»:g 


ARGOl 


deBernardis et al. (1994) 


95 


39.ll 8 8 .7 


ARG02 


Masi et al. (1996) 


95 


46.8^ 


MAX GUM 


Tanaka et al. (1996) 


138 


54.5tf ; 4 9 


MAX ID 


Tanaka et al. (1996) 


138 


46 q+ 21 - 8 


MAX SH 


Tanaka et al. (1996) 


138 




MAX HR 


Tanaka et al. (1996) 


138 


32.7±* 2 9 


MAX PH 


Tanaka et al. (1996) 


138 




Ski 


Netterfield et al. (1997) 


86 


49.0j£° 


Sk2 


Netterfield et al. (1997) 


166 




Sk3 


Netterfield et al. (1997) 


236 


85.0i^ ° 


Sk4 


Netterfield et al. (1997) 


285 


86.0i»;° 


Sk5 


Netterfield et al. (1997) 


348 


69.01™ 


CAT1 


Scott et al. (1996) 


396 


01 -°-13.6 


CAT2 


Scott et al. (1996) 


607 


4Q 1+19.1 



a CMB anisotropy detections reported in publications in 1994 
and later. Many of the new result papers include reanalyses 
of older results. MSAM is not included because of substantial 
spatial and filter function overlap with Saskatoon. 



b Qflat = Qrms-PS \n 3 =l = {j^Y^ 2 STf 
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